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Abstract 

The muon capture rate in the reaction \i~ + 3 He — > Vu + 3 H has been mea- 
sured at PSI using a modular high pressure ionization chamber. The rate corre- 
sponding to statistical hyperfine population of the // 3 He atom is (1496.0±4.0) s _1 . 
This result confirms the PCAC prediction for the pseudoscalar form factors of 
the 3 He - 3 H system and the nucleon. 

PACS: 23.40.-s 



*This paper is part of the thesis works of W. Prymas and N.I. Voropaev 



1 Introduction 



We report on a measurement of the muon- 3 He capture rate to the triton channel, 
\l~ + 3 He — > + 3 H (1.9 MeV), performed to a level of precision unprecedented in 
nuclear muon capture experiments. 

The main information gained from this measurement is the determination of the 
induced nuclear pseudoscalar form factor Fp(q 2 = — 0.954m^) in the A = 3 system 
treated in the frame of the Elementary Particle Model (EPM) pioneered by Kim and 
PrimakofffQ and elaborated recently by Congleton and Fearingp|. The value obtained 
for Fp provides a stringent test of the Partially Conserved Axial Current (PCAC) 
approximation 0, based on the Goldstone mode realization of the approximate chiral 
symmetries of QCD. 

In addition, comparison of the measured capture rate with that calculated according 
to the presently available 0] microscopic model of muon capture by 3 He allows a test 
of the validity of this approach which takes the Meson Exchange Currents (MEC) 
explicitly into account. Using this model, one can evaluate the induced pseudoscalar 
form factor of the nucleon, gp(q 2 = — 0.88m^). Furthermore, our experimental result 
also allows for further tests of the Standard Model in its electroweak sector [§[]. 

In spite of the uncertainties in the various form factors entering the calculation, 
PCAC and EPM predict the capture rate to a precision of 1.4% ||. Until now the 
experimental situation did not equal such a precision, since only three measurements^, 
0, ||, dating back to over thirty years ago were available, with precisions ranging from 
3% to 10%. The situation thus called for a new precision experiment which is described 
in this work. Preliminary results were published earlier 0, 1C . 



2 Experiment 

One of the main advantages in measuring nuclear muon capture on 3 He, compared 
to hydrogen, is the production of a charged particle in the final state, which can be 
detected with high efficiency and good background suppression. For muon capture on 
3 He, the main reaction channel is 

/i-+ 3 He— > 3 H + z/ M (70%), (1) 



accompanied by triton breakup into d + n (20%) and p + 2n (10%)|TT|. 

These reactions occur from the two hyperfine states of the /i 3 He muonic atom, of 
total spin F = and F — 1, see the kinetics scheme in Fig.|l|. Since in the present 
experiment the 3 He target is not polarized and the spin flip is small (cf. results below), 
the hyperfine states are statistically populated, and it is the statistical capture rate to 
the triton channel ^ 

Astat = tAh + tAJj (2) 



which is measured M, 12 



The strategy of the experiment was to select a clean sample of muon stops, iV M , in 
the fiducial volume of an active target well isolated from the target walls and to ensure 
a 100% efficiency in detection of the tritons, N t , from reaction ([!]). Then, the ratio 
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Figure 1: Kinetics scheme of the /i 3 He system. The main disappearance of the muon is 
by decay with the rate Ao (99.7%). The total capture rate from the hyperfine structure 
state F = 0, 1 is Af = A^ + A^, where A^ is the capture rate to the + 3 H channel 
and A^ is the capture rate to the + d + n and + p + 2n channels. 
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Figure 2: Experimental setup and layout of the ionization chamber. 
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Figure 3: (a) Anode layout of the ionization chamber, (b) A typical sequence of anode 
signals registered by the flash ADCs. 



Nt/Np combined with the decay constant of the muon will be a direct measure of the 
muon capture rate. 

The experiment was performed at Paul Scherrer Institute (PSI) using the //E4 
superconducting muon channel. The apparatus is shown in Fig.||]. The target, which 
is also the main detector, consists of a high pressure ionization chamber (IC), filled 
with 120 bar 3 He at room temperature^. The ionization chamber includes a total of 
14 anodes (see Fig.[3|a): 5 entrance anodes (type C), 5 stopping anodes (type B), one 
of which is split, and 3 veto anodes (type A). The anodes have a pad pattern, which 
allows the tracking and location of the stopping muon. The B-anodes cover an area 
of 2.5 x 4.0 cm 2 . The cathode to grid and grid to anode distances are 1.43 cm and 
0.08 cm, respectively. This amounts to a sensitive volume of 15 cm 3 . The cathode 
voltage is — 40 kV, while the grid is kept at — 3.5kV. These conditions result in a 
maximum electron drift time of 3.2 /is. The energy resolution is 30keV (a). The IC 
is surrounded by an array of plastic scintillation counters for the detection of neutrons 
from the breakup reactions and electrons from muon decay. In the present experiment, 
this detector array serves mainly for time scale calibration. A beam telescope of thin 
plastic scintillators provides a fast trigger and also eliminates double-muon events in a 
±8 fis time interval. The muons enter the IC through a Beryllium window of thickness 
4 mm with minimal Coulomb scattering in the window. 

A two-level charge-integrating trigger was used for the B-anodes. The three main 
modes were (threshold 140 keV, typical rate 3000 s _1 ) which fires for each muon 
stopping over the B-anodes, E^ t (two time-separated thresholds at 140 keV, typical 

lr The level of impurities was checked before and after the run using a quadrupole mass spectrometer, 
and found to be 8 x 10~ 5 for 4 He and less than 3 x 10~ 5 for isotopes with Z > 2. 
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rate 10s _1 ) which fires for each muon track followed by a separated triton signal, and 
E t (threshold 1.2 MeV, typical rate lis -1 ) which fires for most tritons and for each 
triton piled up on the muon signal). E^_ t together with E t provide a 100% effective 
trigger for the triton events. For triggered events, the signals from all anodes were 
recorded over a period of 10 /is with flash ADCs (Figfjb). In total, about 4 • 10 8 muons 
entered the fiducial volume triggering E^. The E^ trigger was prescaled by various 
factors k = 500, 1000, and 2000 to reduce the data rate. Finally, 9 • 10 5 prescaled 
and 1.2 • 10 6 E^_ t (+E t ) triggers were registered on tape during a four weeks running 
period. 

3 Data Analysis 

There are two basic event types in this experiment, those in which the triton was 
emitted late enough (At > 500ns) to have its signal cleanly separated from that of the 
muon (see example in Fig.^b) and those in which the two signals overlapped. The bulk 
of the statistics (85%) comes from the separated events. However, to determine the 
muon capture rate, we should measure the total number of tritons including pileups. 
There are two possibilities for measuring the number of the pileups: direct counting of 
the pileup signals or exponential extrapolation to At = of the time distribution of 
the triton signals separated in time, with At = t t — t M , where tt and t M are the arrival 
times of the triton and the muon signals on the IC anodes. 

This data was analyzed using both of these techniques which have different system- 
atic effects and therefore cross-check each other. We label these techniques as analysis 
A (inclusive of pileups) and analysis B (separated events only)|| Furthermore, 
each of these analyses was performed independently in a smaller fiducial volume and 
with a substantially different approach to the systematics: analysis C (inclusive of 
pileups) and analysis D (separated events only)[[y|]. Note that, for instance, analysis 
A and analysis B are identical in all phases except in finding the pileup events. In 
analyses B and D, At was determined with a mid-pulse timing technique, with some 
corrections for the geometrical effects of the muon and triton angles. Analyses C/D 
used the region of muon stops from the center of B2 to the center of B5, while analy- 
ses A/B used the whole region B1-B5. Also, analyses C/D used less vertical space by 
more conservative cuts of the muon drift times. Analysis A, by including pileup events, 
using more B-anodes, and extending the drift volume gave the largest triton statistics 
(1141263 events). 

We should stress that special attention was paid in all of our analyses to guarantee 
strictly identical detection efficiencies for muons with and without tritons. Below, for 
brevity, we present in detail only the results of analysis A together with some remarks 
about and the final results of analyses B/C/D. 

3.1 Selection of Muon Stops 

The position of a muon stop along the beam direction was determined by the energy 
deposits on the stopping anode Bj and the preceeding one. In analyses A/B the muon 
stop position was determined by finding the most downstream anode of those anodes 
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with the biggest number of signals (generally two at most). In analyses C/D this 
position was determined by following a contiguous muon track from a single C-anode 
hit. Vertically, the fiducial volume was determined by cuts on the drift time of muon 
tracks (0.5 [is < t M < 2.3 [is) which set the capture location at least 2 mm from the 
anode and grid, such that the triton track would not touch these electrodes. Note that 
the triton range was 1.52 mm under the conditions of our experiment. Laterally, the 
volume was set by rejecting events in which there was a muon signal on one of the 
A-anodes within ±50 ns from the muon track. This restricted the penetration into the 
A-anodes to less than 0.5 mm. Furthermore, the number of muon stops over the side 
anodes was kept low by triggering only on single, central C-anodes, selecting mostly 
events in the middle of the chamber and further reducing double muon events. 
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Figure 4: Distribution of muon energy depositions S — Ei + 2 • E^i with triggers 
(solid histogram) and E^ t (dashed histogram). The analysis window was set for 
the higher energy peak containing the muons stopping on the anode, while the lower 
energy peak presents muons that passed the anode. 

Finally, we strictly define the characteristics of the muon stop using the muon energy 
deposits Ei on the stopping anode Bi and E^\ on the previous anode. The distribution 
S — Ei + 2 • Ei-i (Fig.p shows a narrow peak for stopping muons. By making a cut 
on S we severely reject any track or noise except a well defined muon stop. This cut 
is loose enough that we estimate from the S'-distribution less than 0.01% difference in 
efficiency for muons followed by capture tritons versus muons without capture reaction. 

3.2 Detection of the Tritons 

In the analysis all triton events were subdivided into the following event types: 
1. Triton signal separated in time from the muon signal. 
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1.1 Triton signal is on the same anode, B i} as the muon stop. 

1.2 Triton signal is split between Bi and an anode neighbour to the one where 
the muon stopped. 

1.3 Triton signal is completely contained on an anode adjacent to the one on 
which the muon stopped. 

2. Triton signal overlaps with the muon signal. 

2.1 Pileup signal is on the same anode, Bi, as the muon stop. 

2.2 Pileup signal is split between Bi and one of the neighbouring B- or A-anodes. 
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Figure 5: (a) Energy spectra of type 1.1 events (triton track contained within a single 
B-anode), (b) of type 1.2 events (track not contained within a single anode) and (c) of 
type 2.1 events (pileups contained within a single anode). Note the vertical logarithmic 
scales. 

Fig.|5| shows the energy spectra of the event types 1.1, 1.2 and 2.1, which dominate 
the overall statistics. The background under the triton peak is due to breakup reac- 
tions. In analysis A the background was subtracted by linear estimates based on the 
background above and below the peaks. The systematic error was then estimated by 
varying the regions of this linear fit. Analyses C/D used polynomial fits resulting in 
comparable errors. The numerical results of analysis A are presented in Table |I|. 

Several corrections were applied to the measured number of N t . The main correc- 
tion (6.445%) is a trivial one taking into account the time beyond the analyzed time 
interval of 5.995 /is. Other corrections compensate mainly for rare event types, e.g. 
tritons passing into one of the A-anodes, not listed in Table 1, and for small losses of 
triton events. The sum of these corrections is 0.735% (see Ref. |10[|), bringing the total 
correction to N t to e t = 7.18%. 




6 



4 Results 



The calculation of the muon capture rate was done according to the formula 

Astat = A • ~ — r r (3) 

k ■ N^(l - e t ) 

where N t is the number of detected tritons, is the number of selected muon stops, 
not followed by muon capture, Ao = 0.45516 ■ 10 6 s _1 is the decay rate of the muon, 
and k is the trigger prescaling factor. Note that in our experiment, the muon stops 
which are not followed by muon capture are the events triggered by only. This 
trigger contains practically none of the muon capture events (except for about 10% of 
low energetic breakup events giving a negligible correction to iV M of 1.4- 10~ 4 ). 
The statistical error is determined by the following expression 




^Astat - ^stat ■ 1 77 H 77 (4) 



where F is the fraction of muon stops used in the analysis. (In the approximation of 
Eq. [| accounting for all muons would leave only the statistical error of N t .) In the 
considered case of analysis A: N t = 1141263, = 349479, F = 0.356, while the total 
number of selected muon stops corrected by the prescaling factor is k-N^ = 374028500. 
Substituting these numbers for each k value separately in (|3|) and (|j) and combining 
results statistically we obtain 

A s tat = (1496.2 ±2.6 stat )s- 1 . (5) 

The systematic error was evaluated by varying the energy and time cuts in wide 
ranges and comparing results from independent data on different anodes and under 
different running conditions. The main contribution to the total systematic error comes 
from uncertainties in the background subtraction (0.16%). The residual systematic 
errors in the determination of iV M and Nt are 0.03% and 0.05% respectively, giving a 
total systematic error of 0.17%. 



Event type 


N t 


N t /Nl ot (%) 


BG 


BG/N t (%) 


1.1 


843322 


73.9 


30882 


3.7 


1.2 


77881 


6.8 


11943 


15.3 


1.3 


52638 


4.6 


2684 


5.1 


2.1 


157165 


13.8 


20144 


12.8 


2.2 


10257 


0.9 


9889 


96.4 


sum 


1141263 


100.0 


75542 


6.62 



Table 1: Total statistics of analysis A. Nt is the number of triton events found on 
each anode, N^ ot the sum, and BG is the number of evaluated background counts 
underneath the triton line. 
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Table || presents the partial results in measuring A sta t (analysis A). These results 
show a remarkable stability. The final result from analysis A is 

analysis A : A stat = (1496.2 ± 2.6 stat ± 2.6 syst ) S - 1 (6) 

The results from the other three analyses are: 

analysis B : A stat = (1495.8 ± 3.0 stat ± 2.8 syst )s" 1 (7) 

analysis C : A stat = (1497.4 ± 3.1 stat ± S^^s" 1 (8) 

analysis D : A stat = (1496.0 ± 4.5 stat ± 3.0 syst ) S - 1 (9) 

Let us recall, that the principal difference between analyses B/D and A/C is that 
B/D used only the separated events, requiring extrapolation back to At = 0. This 
method relies on the precision of the determination of the absolute value of the At = 
position. We have shown that the mid-pulse to mid-pulse technique determines this 
position with a precision better than ±2 ns. In analysis B this was proven using the 
neutron detectors to compare the arrival times of decay electrons with neutrons from 
the breakup channels and correcting for the estimated time of flight. In analysis D, 
At = was confirmed by checking the average angle and curvature of muon tracks 
to estimate any offset to the mid-pulse technique. In all analyses, several other sys- 
tematic effects were checked. These include: edge effects, binning errors, triton track 
angle effects, and pedestal shift effects. The very good agreement among all analyses 
demonstrates the low level of remaining systematic errors. 

Finally, we present our result as 

A 9 t a t = (1496.0 ± 4.0)s- 1 (10) 

where the error is the combined statistical and systematic error of the larger statistics 
analyses A and B. 

The measured value of A sta t corresponds to a statistical population of the singlet 
and triplet states of /x 3 He atoms. An assumed hyperfine transition from the singlet 
to the energetically lower triplet state with a rate A ^ 1 = 0.01 ns~ l would change the 



k 


2000 


2000 


1000 


500 




53664 


70619 


25729 


199467 


N t 


327252 


429989 


78622 


305400 


Bl 


1497±16 


1500±15 


1514±23 


1509±11 


B2 


1483±15 


1489±14 


1494±21 


1482±10 


B3 


1484±15 


1471±13 


1482±20 


1503±10 


B4 


1498±14 


1486±13 


1514±20 


1516±10 


B5 


1506±14 


1525±13 


1497±20 


1492±10 


B1-B5 


1493±6 


1495±5 


1501±8 


1500±4 



Table 2: k, N^, N t and partial results on A sta t from different B-anodes and run groups 
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Figure 6: Time distribution of separated fi—t events. The events near t = are absent 
due to muon-triton signal overlap. In the time region for analysis B (0.6 — 6.05 /xs) a 
pure exponential slope with A = (0.4573 ± 0.0005) s _1 was found, in accordance with 
the known lifetime of the /x 3 He muonic atom. 



observed A s t a t by +0.17%. Fortunately, we can limit this process in our experiment by 
comparing the observed time distribution (Fig. 6) of the tritons with the expected one: 



dNt 
dt 



-(Ao+AO+A - 1 )* 



-A H e 



3 U 



+ 



A 



0^1 



Ah 



4 A - 1 + A° - XI 



+ 

(Ao+AO+A - 1 )^ 



(11) 



Using the rate A sta t = (A H + 3A H )/4 determined above and the total spin averaged 
capture rate (A° + 3Aj)/4 = (2216 ± 70) s -1 |T4| together with the theoretical values 
for the differences between the singlet and triplet capture rates A^ — A^ = 578 s _1 |f2[ 
and A° — Aj = 450 s" 1 [JTS]] we obtain 



A 



0^1 



(0.006 ± 0.008) fis- 1 



(12) 



thus restricting the possible influence of this process to a level of <0.14%. 

Another uncertainty might arise if a significant fraction of the 2s-state were metastable 
with a lifetime comparable to the muon decay rate. The total initial population for 
the 2s-state is 7%[0J. According to the present understanding|fr6|, [T7|, however, the 
quenching of the 2s-state occurs on a much shorter time scale and does not influence 
our result. But again, we can limit this effect with our own data by comparing the 
number of tritons, N t , obtained by the extrapolation method (analyses B/D) with that 
obtained by direct counting of pileups (analyses A/C). From this comparison we obtain 
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the upper limit for the lifetime of the 2s-state 

r 2s < 50 ns. (13) 
The effect of this transition is indeed negligible. 

5 Discussion 

The obtained result can be analyzed in the framework of the EPM0. In this model, the 
weak current of the 3 He- 3 H transition is parametrized by four form factors Fy, Fu, F\, 
and Fp, evaluated at the relevant value of the four-momentum transfer q$ = — 0.954m^, 
see e.g. Ref.|§. Two of these parameters, the vector and magnetic form factors, Fy 
and Fm, are derived from the conserved vector current (CVC) theorem and from the 
results of electron elastic scattering by 3 He and 3 H[0: 

F v (q^) = 0.834 ± 0.011, F M (q%) = -13.969 ± 0.052. (14) 

Taking these data as inputs, our measurement determines the correlation between the 
pseudoscalar form factor Fp(gg) and the axial form factor Fa(<2o), as shown in Fig. 7. 
The width of the allowed region is determined mostly by the error in Fy(q^). In a 
further step we use F A (q$) = 1.052 ±0.010 extrapolatedg from F A (0) = 1.212 ±0.005 
as determined from the 3 H beta decay rate. Note that this extrapolation involves 
some theoretical consideration based on the impulse approximation that results in the 
increase of the error in Fa(c/o). This procedure allows determination of the pseudoscalar 
form factor 

Fp(ql) = 20.8 ± 2.8. (15) 

The main contribution to the error comes from the error on F^q^), while the contri- 
bution from our measurement is only ± 0.5 . 

The obtained value of Fp can be compared with the PCAC prediction that relates 
the pseudoscalar and axial form factors: 

PCAC = m M (M3 He ±M3 H ) 2 
m z _ gz 

This relation is shown in Fig.|7j by the dot-dashed line. Taking again the extrapolated 
value for Fa, one obtains 

^p CAC (9o) = 20.7 ±0.2 (17) 

in perfect agreement with our experimental result (15). 

The availability of precise three-body wave functions of the A = 3 system makes 
it possible to perform calculations of the muon capture rate on a "microscopic" level. 
Calculations done in the impulse approximation lead to a capture rate0] of A sta t = 
1304 s -1 which deviates from the observed rate by 13%. This clearly indicates that 
MEC must be accounted for. Full microscopic calculations including MEC have been 
performed recently by Congleton and TruhhkQ. The authors obtained A stat = 1502 ± 
32 s -1 in excellent agreement with our experimental result (10). This agreement can 
be interpreted as a confirmation of the validity of the approach used to include the 
MEC in the microscopic theory. 
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Figure 7: Plot of axial form factor Fa versus pseudoscalar form factor Fp showing the 
constraints by the experimental errors of our measurement (full lines) and with errors 
in F v and F M added (dashed lines). The vertical dotted lines give the constraints on 
Fa, and the dot-dashed line presents the PCAC prediction. 

Note that these calculations required the knowledge of the nucleon form factors gy, 
<?A, <?m, and gp at q%. Three of these, namely gy, gu, and g&, are experimental values, 
while the pseudoscalar form factor g P was determined by the PCAC relation || to be 

flP^fog) = 8.12 . (18) 



In this context, we wish to point out two recent publications fll8|, where, for the 



nucleon, a formula identical to the PCAC relation was derived within the framework 
of chiral perturbation theory. Assuming the validity of the microscopic model within 
the errors cited in Ref. |§], we can determine gp(<2o) from comparison of the calculated 
value of A s tat with our experimental result: 

g P (q 2 ) = 8.53 ±1.54. (19) 

The error of this result is largely dominated by the errors in the calculation of the 
theoretical muon capture rate and not in the precision of the present experimental 
result. Note also that our value of <?p(<Zo) agrees well with the PCAC prediction 



1.05 ±0.19. (20) 



while the recent experiment |20| on radiative muon capture on hydrogen shows a sig- 
nificant deviation. Further studies of muon capture on hydrogen are definitely needed 
for a determination of g-p with higher accuracy. 
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